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ABSTRACT 


Various applications of attitude estimation as applied to the opera- 
tion of the Orbiting Astronomical Observatory are enumerated and ex- 
plained in some detail. The techniques used, the results of these tech- 
niques, and the problems encountered during the first nine months of the 
mission are delineated. The operation of the Orbiting Astronomical 
Observatory including the spacecraft and the supporting ground system 
is briefly described. 


INTRODUCTION 

On December 7, 1968 the NASA's Orbiting Astro- 
nomical Observatory (OAO) was placed into a near cir- 
cular orbit inclined 35 degrees from the equator by an 
Atlas-Centaur booster. The OAO, orbiting 480 miles 
above the earth's surface, doubled existing stellar data 
clear of the earth's absorbing atmosphere in its first 
thirty days of operation. After nine months of suc- 
cessful operations this first astronomical observatory 
in space has made over 3000 observations of stellar 
objects and has attained some 5000 different attitudes 
distributed over the entire celestial sphere. 

Two experiments are aboard the OAO spacecraft: 
One designed by the University of Wisconsin and the 
other by the Smithsonian Astrophysical Observatory 
(SAO). The Wisconsin Experimental Package (WEP) 
includes four stellar photometers, two scanning spec- 
trometers, and one nebula photometer. The primary 
function of the WEP is to gather spectral energy dis- 
tributions on selected stars and nebulae in the ultra- 
violet range of 1000 to 4000 Angstrom (A). The SAO 
experiment contains four 12-inch telescopes each of 
which images a star field onto an ultraviolet sensitive 
photocathode. The scientific objective of SAO is to 
measure the brightness of many celestial bodies in four 
separate ultraviolet spectral bands ranging between 
1100 A and 3000 A. 

THE SPACECRAFT AND ITS CONTROL 

The OAO spacecraft is an octagonally shaped 
aluminum structure ten feet long and six and two-thirds 


feet across any two parallel sides of the octagon. There 
is a hollow central tubular area four feet in diameter 
running the length of the spacecraft which is used to 
house the experimental equipment. Two hinged sun- 
shades approximately four feet square are attached to 
each end of the spacecraft to protect the experimental 
equipment from sunlight. Attached to each of two op- 
posite sides of the octagon at an angle of 34 degrees 
from the longitudinal (optical) axis are four solar cell 
paddles whose combined dimensions are nearly eleven 
feet by nine feet. Solar cells are mounted on both sides 
of all eight paddles. The major subsystems of the 
OAO are stabilization and control, data processing, 
communications, and power. Thermal control is pas- 
sive. The entire spacecraft weighs over two tons 
which includes 1000 pounds of experimental equipment. 

The primary mode of control is by stellar guidance 
so as to point either of the experiments to within one 
arc-minute of a specified target star. Control perpen- 
dicular to the optical axis is also required for maneu- 
vering, power, thermal, and shading considerations. 

This three dimensional control is accomplished by us- 
ing various combinations of six orthogonally mounted 
gimballed startrackers each having a field of view of 
one degree and an excursion of approximately forty 
degrees in any direction. A seventh tracker, bore- 
sighted with the optical axis, is also available which 
provides single axis (optical) control. Averaged error 
signals from the trackers drive three orthogonally 
mounted "fine” momentum wheels which absorb the 
extraneous torques that would cause attitude errors. 

The momentum in the fine wheels iB continuously 
"dumped" by a magnetic unloading system consisting 
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of three flux gate magnetometers, three torquer coils 
(electromagnets), and a processor. The momentum in 
the fine wheels may also be removed by a gas jet system. 

The OAO may also be controlled by three rate 
gyros which actuate either the fine momentum wheels 
(Hold On Wheels) or the gas system (Hold On Jets). 

The maximum drift rate for any one gyro as obtained 
from in-flight data is about 0.2 degrees per hour. 

The spacecraft is maneuvered from one attitude 
(two axes) to another by means of three "coarse" iner- 
tia wheels in an open loop manner. Two or more 
trackers are required to track continuously during the 
slew and averaged errors from preset gimbal angles 
settle the spacecraft when the slew has terminated. 
Slewing plus settling times vary with the axis and angle 
of rotation but is usually better than three degrees per 
minute. Small closed loop slews (several degrees) may 
also be accomplished with the fine wheels by changing 
the commanded gimbal angles. 

The OAO has two special modes of control: "sun- 
bathing" and "sun-pointing". In the former, sun sen- 
sors are used to actuate gas jets which align the nor- 
mal of the solar cell array to the sun when in sunlight. 
Rate gyros then hold this position while in darkness. 
Sun-pointing is similar except in this case the optical 
axis is pointed toward the anti- sun line. This is the 
pointing the spacecraft automatically obtains upon sep- 
aration from the launch vehicle. 

In addition to the sensors mentioned above there 
are four sun sensors mounted so a6 to give solar as- 
pect data at any attitude. It is the purpose of these 
sensors to provide independent attitude information 
while in sunlight. 

The data processing subsystem can accept 168 dif- 
ferent control commands in addition to the gimbal and 
slewing commands. All commands may be executed in 
real time or stored (up to 256 commands) and executed 
at a specific time later as a function of an on-board 
clock. 

THE GROUND SYSTEM 

For the most part, the actual functions performed 
by the OAO are initiated by one of the experimenters 
although several major elements of the ground system 
lie between his request and the commands executed by 
the spacecraft. The experimenter specifies a sequence 
of target stars (pointings) and the associated experi- 
mental equipment commands necessary to obtain the 
desired scientific objectives. This information (Experi- 
menter's Target List) is passed on to the Mission 
Computing Group where it and other inputs are entered 
into a complex computing system known as the Support 
Computer Program System (SCPS). The SCPS (which 
resides in a large scale-high speed computer) deter- 
mines the total attitude as a function of the target and 


other geometric considerations. From this, the num- 
ber of turns for each inertia wheel required to slew 
from the previous attitude to the present one is then 
determined. The gimbal angles for each tracker and 
their on-off schedule as a function of occultations by 
sun, moon, or earth are also generated. One of the 
final outputs of the SCPS is the ordered list of com- 
mands in spacecraft format which are to be executed 
and their time of execution. Based upon the schedule 
of ground contacts with the spacecraft, the number of 
commands to be executed, and the current memory as- 
signments the SCPS also determines when and where 
the commands are to be loaded into the spacecraft 
memory. 

The image of the OAO command memory along 
with certain ground procedures generated by the SCPS 
is then routed to the OAO Operations Control Center 
(OCC) which in turn transmits this "contact message" 
to one of five remote sites. These sites are: Rosman, 
North Carolina; Quito, Ecuador; Santiago, Chile; 
Tananarive, Madagascar; and Orroral, Australia. At 
the scheduled contact the remote station establishes 
communications with the OAO, loads new commands in 
memory (if necessary) , and gathers telemetered data 
(real time and stored). The OCC monitors these con- 
tacts (in the case of a Rosman contact it also replaces 
the remote station functions) , and displays the return- 
ing telemetered data. During the contact, real time 
commands already at the remote site may be executed 
and real time commands from Goddard Space Flight 
Center (where the SCPS and OCC are located) can be 
received by the site and relayed to the spacecraft, 

ATTITUDE ESTIMATION APPLICATIONS 

The applications for attitude estimation techniques 
in the OAO program range from coarse estimates 
(several degrees) from the solar and magnetic sensors 
to a precise attitude determination (seconds of arc) 
from as many as six startrackers. 

A coarse attitude estimate is used as an inter- 
mediate step in obtaining stellar guidance. The transi- 
tion between the coarse attitude estimate and stellar 
control is obtained by generating a sequence of star- 
tracker gimbal angles consistent with a sequence of 
attitudes. This sequence of attitudes sweeps out the 
most probable region associated with the coarse atti- 
tude estimate. Stellar guidance is then obtained when 
two or more trackers acquire their preassigned stars 
during the search maneuver. 

The coarse estimate of the OAO's attitude is ob- 
tained by placing the spacecraft in an attitude hold using 
the rate gyros and fine momentum wheels. Sun sensor 
and magnetometer data are then collected while the 
spacecraft is in contact with a ground station (approxi- 
mately ten minutes). This data is then relayed, by high 
speed data links, to the control center where it is 
passed to the SCPS. Since the attitude of the spacecraft 
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is fairly stable during the collection of the data (rela- 
tive to the accuracy of the data), it is assumed that all 
the data refers to the same attitude and is thus used 
collectively to obtain a weighted least squares estimate 
of th t attitude. The star search commands, based upon 
the attitude estimate, are then generated by the SCPS 
and transmitted (via the total ground system) to the 
spacecraft at some subsequent contact in real time. 

Even though the magnetometers were not calibrated 
and give noisy data (their original purjjose was solely 
to unload momentum from the fine wheels) the above 
procedure yields an adequate solution and has been used 
to place the OAO under stellar control approximately 
thirty times in a nine month period. The success of 
this technique has saved valuable time and extended the 
life of the OAO's gas supply since it is applicable at 
any orientation whereas the original design concept re- 
quired a reorientation to align a particular axis with 
the sun. Stellar control would then be accomplished by 
rolling about this axis until two or more startrackers 
simultaneously acquired stars. 

The OAO's capability to hold an unknown attitude 
under gyro control has also been used to derive an es- 
timate of attitude using only magnetometer data. To 
date such rough estimates have had very limited appli- 
cation. However, on one occasion such an estimate was 
used to verify that the telemetry from a sun sensor had 
been misinterpreted. With calibrated magnetometers 
this capability could provide useful attitude information 
while the spacecraft is in darkness. 

When under stellar control, attitude determination 
techniques using startracker data have been used for 
an independent verification of the commanded attitudes 
and to determine the accuracy of the various sensors. 
Least squares solutions using all tracker data (usually 
from two to four trackers) as an aggregate and in vari- 
ous combinations provides an excellent means for 
evaluating the performance of each tracker and whether 
any misalignments exist between them. 

More sophisticated techniques have been employed 
to determine the values of the tracker misalignments 
from inflight data. This has been necessary since 
shifts due to launch stresses and other effects have 
been greater than expected. By use of this procedure 
the pointing accuracy has been better than the design 
specification of one arc-minute. 

Attitude determination from tracker gimbal angles 
is also performed in the processing of the experimental 
data. This allows a correlation of the data to a known 
star and aids in the evaluation of the scientific results. 

Stellar control lias provided the opportunity to study 
the correlation between the errors of the sensors and 
the resulting error in attitude. In this manner a good 
deal of experience has been obtained in the selection of 
weighting factors for the various sensors as well as 


interpreting the results based on residual errors and 
the angular separation between the solar and magnetic 
vectors. A technique has recently been developed which 
will determine from in-flight data any systematic errors 
caused by: the misalignment of the magnetometers, 
the permanent magnetic moment of the spacecraft, and 
the effects of the torquer coils on the magnetometers. 

It is felt that this information will help eliminate the 
present wide dispersion in the magnetic data. 

THE ATTITUDE ESTIMATE 

All of the attitude estimates referred to in the 
preceding sections are obtained from the same mathe- 
matical algorithm which yields a weighted least squares 
estimate of attitude. The only difference between the 
various applications is the selection of the weights 
and/or the type of data to be used. 

The method (delineated in an appendix) is based 
upon two different vector parametrizations of three- 
dimensional rotations. One vector has been referred 
to as the Gibbs vector and the other's components are 
three of Euler's symmetrical parameters. The direc- 
tion of both vectors defines the axis of rotation and 
their lengths are trigonometric functions of the angle 
of rotation. One vector has a length equal to the tan- 
gent of the half-angle of rotation and is denoted as the 
Y vector. The other vector is denoted as Z and has 
length equal to the sine of the half-angle of rotation. 

The Z vector is used to obtain a vector expression 
for the "smallest" rotation which will align an estimate 
of a rotated vector to its true value. A generalized 
weighted least squares criterion is thus established by 
requiring the sum of squared lengths of all such Z's 
(premultiplied by a symmetric weight matrix) to be a 
minimum. The resulting equation is then simplified 
by assuming that the weight matrix is the identity times 
a scalar (this implies that the component errors of 
each Z vector are equal and independent). 

The Y vector is used to express the least squares 
condition in vector notation where the only variables 
are the three components of the Y vector. The neces- 
sary conditions for an extremum are then applied to 
this function and the result is a vector equation in terms 
of the Y vector. The least squares solution can then be 
obtained by finding the largest zero of a fourth degree 
polynominal and then solving a linear system of three 
equations. In practice, this approach is not taken since 
the vector equation is also amenable to a simple suc- 
cessive substitutions iteration which converges rapidly. 
Thus an approximate solution is constructed to start 
the iteration. 

When only two data vectors are given the least 
squares rotation can be expressed as a product of two 
rotations, each of which are obtained by juxtaposing 
vectors that are linear combinations of the given vec- 
tors. The scalar coefficients of these linear combinations 
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are explicitly given as functions of the lengths and dot 
products of the given vectors. This method yields the 
attitude estimate when only two data vectors are avail- 
able and provides the first approximation when more 
than two vectors have been measured. 

RESULTS AND PROBLEMS 

Attitude error signals from the experiments are 
not available on the present OAO. Therefore, it is dif- 
ficult to ascertain the absolute error in attitude while 
under stellar control. However, attitude determinations 
using telemetered gimbal angles from several trackers 
over a ten minute interval usually yields residual er- 
rors of one arc- minute or less (the smallest command 
increment is twenty arc-seconds). To maintain this 
precision, however, it has been necessary to occas- 
ionally re-evaluate the misalignments between the 
trackers. More elaborate misalignment models are 
being considered to eliminate this miror problem and 


to improve the overall accuracy. Since the misalign- 
ments have been as good as thirty arc-seconds over ex- 
tended periods it is felt that much of the variation can 
be accounted for with additional parameters in the 
model. 

In the following table are the results of fifteen 
typical attitude estimates using only sim sensor and 
magnetometer data. Column one contains the number 
of data points used in the estimate-each data point gives 
two vectors (solar and magnetic). The second column 
shows the angular separation between the solar and 
magnetic vectors. The length of time over which the 
data were taken is given in column three. Columns 
four and five give the actual angular error of the solar 
and magnetic data respectively (the attitude obtained 
from tracker data serves as the reference). The next 
two columns displays the solar and magnetic errors 
using the solar- magnetic estimate as a reference. The 
last column contains the angle of that rotation necessary 


No. of 
Data Points 

Angular 

Separation 

Duration 
min: sec. 

Actual 
Sim Error 

Mag. Error 
Mean, Max. 

Residual 
Sun Error 

Mag. Residual 
Mean, Max. 

Attitude 

Error 

G 

36 to 38 

0:47 

0.02 

3.9, 4.9 

0.18 

3.9, 4.8 

1.1 

G 

62 

64 

1:03 

0.25 

6.4, 8.6 

0.32 

5.8, 8.0 

1.1 

10 

39 

44 

1:50 

0.14 

3.6, 4.5 

0.20 

3.4, 4.4 

1.2 

11 

37 

29 

9:42 

0.06 

2.6, 4.0 

0.08 

2.2, 3.3 

1.8 

11 

56 

49 

2:21 

0.06 

4.7, 6.0 

0.37 

4.4, 5.9 

1.8 

14 

81 

82 

4:11 

0.08 

7.1, 9.2 

0.31 

6.8, 9.9 

2.8 

16 

40 

50 

3:25 

0.26 

2.2, 3.5 

0.15 

2.3, 3.8 

0.9 

17 

37 

28 

3:56 

0.06 

3.4, 5.4 

0.16 

2.6, 5.3 

3.5 

18 

47 

25 

8:23 

0.02 

1.4, 3.3 

0.05 

1.4, 3.3 

0.5 

24 

69 

45 

7:52 

0.06 

4.5, 6.3 

0.30 

4.2, 6.0 

0.3 

25 

80 

107 

7:07 

0.35 

5.4, 6.4 

0.31 

5.4, 6.4 

0.8 

25 

88 

78 

7:52 

0.09 

6.1, 8.1 

0.32 

5.8, 8.1 

1.0 

25 

29 

41 

6:49 

0.22 

1.7, 3.4 

0.08 

2.2, 4.0 

1.5 

25 

61 

46 

6:34 

0.21 

6.1, 8.5 

0.33 

5.4, 8.4 

1.9 

25 

74 

61 

6:49 

0.26 

5.5, 8.1 

0.30 

4.7, 7.6 

2.1 


All errors and angular separation are in degrees. 

In all attitude estimates the sun sensor data was given a weight of 4 and the magnetometer data a weight of 1. 
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to bring the solar-magnetic estimate in agreement with 
the stellar determination. These figures clearly show 
thr, dispersion in the magnetic data that was referred 
to previously. The errors range from as good as our 
capability to predict the magnetic vector (the time of 
the data is only known to within 1G seconds) to a factor 
of five times worse. This makes it very difficult to 
determine the accuracy of the estimate in practice, 
especially when the angular separation is poor (less 
than twenty degrees). In this case, small residual er- 
rors do not necessarily mean a good estimate. On the 
other hand, a separation of less than twenty degrees 
may yield a good estimate if the magnetic data are 
good. 

It is rather apparent from the table that attitude 
estimates based upon the magnetic data only would be 
quite poor in general. However, we made such an es- 
timate for what was considered to be the best case (the 
one with eighteen data points) merely to see what type 
of accuracy could be obtained by this method. The at- 
titude error in this case was 0.95 degrees. Thus, if 
the larger errors in the magnetic data can life elimi- 
nated, useful attitude estimates could be obtained while 
the spacecraft is in darkness. 

Of the thirty some star searches attempted, only 
three have failed due to poor attitude estimates. The 
present search capability', due to a ground constraint, 
is about 3.5 degrees which is the maximum error of the 
fifteen cases discussed earlier. However, in the fail- 
ure cases data was limited to a small interval of time 
where the angular separation was also small. Two 
searches have also failed because one or more track- 
ers locked onto wrong stars (a bright neighbor of the 
selected star). Normally, logic aboard the OAO pre- 
vents this by requiring simultaneous acquisition of 
stars. This logic which is set whenever all trackers 
lose their stars, has been by-passed several times by 
ground commands when it was apparent that stellar 
control had been lost (only one tracker tracking). This 
action is justifiable in that it also prevents further loss 
of control. 

The occasional large errors in the magnetic data 
have alerted us to a potential problem in estimating 
attitude from multiple data representing essentially 
only two directions. When one direction is weighed 
heavily over the other the solution is weak in one di- 
mension geometrically and in another statistically. The 
net result is that there is another solution (a saddle 
point) to the necessary conditions for extremum near 
the desired least squares solution. This is especially 
troublesome for small angular separations. 

CONCLUSION 

In every respect, the present OAO mission must be 
considered a resounding success. Although the OAO is 
the most complex scientific satellite ever built it has 
operated almost flawlessly for nine continuous months. 


The ground system has been able to keep the satellite 
productive around the clock and both experimenters 
are continuing to receive valuable scientific data. The 
problems, thus far, have been minor and solved with 
work-around procedures on the ground. 

With future OAO spacecraft carrying such equip- 
ment as integrating as well as rate gyros, ultra pre- 
cision experiments, and an on-board computer it ap- 
pears that the GAO program will continue to offer 
challenging and exciting possibilities in the field of at- 
titude determination and control. 
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APPENDIX 


A VECTOR APPROACH TO WEIGHTED LEAST SQUARES ATTITUDE ESTIMATION 


Given the angle of a rotation, , and the rotation 
axis defined by the unit vector X it is well known that 
the matrix, R, of the rotation can be expressed as 


Z 


/ « I . 

A (1+U n T n) " " 


R x (i9) - cos 6 I ♦ (1 - cos tf) XX T - sin e» X - (1) 

Here X denotes the skew-symmetric matrix formed 
from the components of X such that for any vector V, 

XV r X » V (I is the 3 *3 identity matrix and the super- 
script T denotes transpose). Furthermore, any ortho- 
gonal matrix with determinant equal to plus one can be 
expressed in the form of Eq. (1). Except when other- 
wise stated, it is assumed that X has been selected so 
0 < 5 < 7 7 . 

If we define the vectors Y and Z as 


where U = |U| U n and T - Tl T n . In other words, the 
angle of the rotation is the angle between U and T, 
and the rotation axis is perpendicular to the plane con- 
taining these two vectors. If |U| T| then U = RT. 
In particular, if U is a measurement (including error) 
of vector components relative to a local right-handed 
orthonormal coordinate Bystem, T is the same physical 
vector but relative to a fixed reference system (also 
right-handed and orthonormal), and R is an estimate 
of the rotation (attitude) relating the two systems then 
the vector error that can be corrected by a r otation is 
given by 


Y = t«n(0/2)X, Z = sin(6'2)X. 


AZ = 


^2(1 +U n • RT n ) L '" ' RT " ' 


( 2 ) 


then we find, by standard trigonometric identities, that The quantity U - RT is the error in a true vector sense— 

R may be written as it cannot be removed, in general, by a rotation. 

For each such measurement, U t , there corresponds 
a AZj. Thus, the weighted sum of "rotational errors" 
squared is given by 


R = 


1 * Y • Y 


[(1 - Y • Y) I + 2YY T - 2y] 


or 

R = (1 -2Z • Z) I + 2ZZ t - 2 /T^T^Z Z . 

In this manner, the vector components relative to a 
rotated system can be expressed, by vector relations, 
as a function of Y or Z and the components relative to 
a fixed system. For every rotation there exists a Z 
vector, uniquely except for 180 degree rotations, A 
unique Y vector exists for any rotational matrix whose 
trace is different from minus one. 

With the above definitions it then follows that if 
T and U are any two vectors, the rotation, R, which 
will align U and RT with the ir.iaunum angle of rotation 
can be expressed by the Z rotation vector as 


f(R) = 2^ ( p . Az .) Tp . AZ i • (3) 

i* 1 

where m is the number of measured vectors and P, is 
the weight matrix for the i ,h measurement. If P, r p, I 
(p j a scalar) and equation (2) is substituted into equa- 
tion (3) we obtain 


f(R) = ^ ( W t ” RV l)* ’ (4) 

i* I 
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where 


and 



_ U J_ 

lu.l 


The matrix A T A is symmetric and positive- semi- 
definite. Let d* dcnotcB the i ,h ordered (dj 2 > Hj 2 
> tl > 0) non-negative eigen.alue of A T A correspond- 
ing to the normalized eigenvector N,. The vectors 
AN, ( i = 1,2,3) then constitute an orthogonal system 
ofvectors with I AN, d,. Thus, an orthonormal sys- 
tem of vectors T, ( i = 1,2,3) can l>c constructed 
such that AN, =d, T,, If d, " 0 then T, may l>e any 
arbitrary unit vector; otherwise 





T , 



The function f(R) , the weighted squared errors, 
is minimized when the function 


g(R) = • RV, 

i* i 


Tj and T 3 are then similarly constructed. Let N and 
T denote the orthogonal matrices obtained by juxtaposing 
the vectors N, and T, respectively. Then from the 
constructions above it follows that T* 1 AN ' D , where 
D is a diagonal matrix with entries d, (ordered). In- 
troducing the orthogonal matrix M = N* 1 RT ( R is the 
desired least squares rotation) we find th-.t the expres- 
sion tr(KA) as a function of M becomes 


is a maximum. This expression may be written as 


g(R) = (RV,) W, T 

»« 1 


2 ^ tr(RV, W, T ) 
»■ 1 


tr *22 V ,*- T 
_ 1*1 


tr(RA) 


tr(RA) = tr(NMf A) 


tr(MD) 



!• l 


For M orthogonal and d, >0, the maximum value of 
tr(RA) is obtained when M = I. Since M N R T , 

R will be a rotatior matrix when M r I if ! N | |T| - l. 

If |N T = -1 then the maximum value of tr(R A) for 
iRl = 1 is obtained when 


where tr denotes "trace of and 



i • t 


Hence, the least squares rotation is obtained by maxi- 
mizing the function g(R) r tr(RA). 

Except for notation, Eq. 4 is the same as the equa- 
tion in Ref. 4 to be minimized and the solutions there 
may also be applied to Eq. 4. There is, however, con- 
siderable difference in the definitions of the vectors 
used in the equations. As applied to attitude estimation, 
the vectors in Ref. 4 are direction cosines (although the 
solution also applies to unnormalized vectors) whereas 
here they are generated from the measurements (nor- 
malized or not) so a6 to minimize the weighted rotation 
error. Tne following solution to minimizing Eq. 4 is 
similar to that given in Ref. 4. 


M 



0 

1 

0 



This latter case is equivalent to changing the sign of 
Tj in the definition of T. Thus, R = NT 0 t , where 
T 0 = (Tj.Tj. N T Tj ) is the desired solution to the 
least squares condition. This solution can easily be 
generalized for vectors of arbitrary dimensions. 


When only two measured vectors are given it can 
be shown that 




V), + xWj 


W, ♦ x»j 


xV - V 

l*v, -v,| 


* V i 4 V z 

iyv, ♦ v a l 
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N a - N, « Nj , and T 3 T, - T 2 gives the least squares 
solution if * and y satisfy the condition x 4 y > 0 as 
well as the simultaneous equations: 


(V, ' Vj) xy 4 (V | ' Vj) ( x - y ) - V, • V, 0 
(*, ' *,) 4 (W, • * 2 ) fx- y) - xy Wj • W 2 = 0 
The solution to this pair of equations is given by 


If Eq. 4 is written as a function of the Y rotation 
vector and the conditions for an extremum applied, 
the resulting three scalar equations can Ire expressed 
by a single vector equation as 


• v .) • Y< ( v . • y )(*. • Y ) • », 


n i /a 2 4 4bc 
x * 2c • 

■at 4 4bc 

y = Jc • 

with 

» = ( v j • v 2 ) (v f a) - ( v , • v .) (V*,) • 

b - (v,-v 2 ) (*,•*,) 4 (»,•*,) (v a -v 2 ) . 

c = (w, • w 2 (V, • V 2 ) 4 (V, • v 2 ) (w 2 • w 2 ) . 


<I 4 Y.Y,£[(V, * Y ) w , 4 (*1 ■ Y) v, 4 W, X V,] . 


Ref. 2 gives a solution of this equation which requires 
obtaining the largest zero of a fourth order polynomial 
and then solving a linear system of three equations. It 
has also been solved successfully by a successive sub- 
stitutions iteration (divide both sides of the equation by 
the coefficient of Y on the left side). In this case, the 
V ( are first rotated by an approximation of R. The 
rotation corresponding to Y is then the correction to 
the approximation necessary to obtain the least squares 
solution. 
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